Dipolar Bose-Einstein condensate in a ring or in a shell 
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We study properties of a trapped dipolar Bose-Einstein condensate (BEC) in a circular ring or a 
spherical shell using the mean-field Gross-Pitaevskii equation. In the case of the ring-shaped trap we 
consider different orientations of the ring with respect to the polarization direction of the dipoles. In 
the presence of long-range anisotropic dipolar and short-range contact interactions, the anisotropic 
density distribution of the dipolar BEC in both traps is discussed in detail. The stability condition 
of the dipolar BEC in both traps is illustrated in phase plot of dipolar and contact interactions. We 
also study and discuss the properties of a vortex dipolar BEC in these traps. 

PACS numbers: 03.75.Hh,03.75.Nt 



I. INTRODUCTION 

After the experimental realization of dilute trapped 
Bose-Einstein condensate (BEC) of alkali-metal atoms 
PP, it was realized that for attractive atomic contact in- 
teraction the BEC is stable for the atomic interaction 
(measured by the scattering length) and number of atoms 
below a certain limit depending on the trap p~H3] . In the 
case of 7 Li atoms with attractive atomic interaction, in 
the harmonic trap used in experiment [2 j, the conden- 
sate was stable for less than about 1400 atoms. For re- 
pulsive contact interaction, the trapped BEC is uncondi- 
tionally stable for all values scattering length and number 
of atoms p. 

More recently, there has been experimental observa- 
tion of BECs of 52 Cr pQ [5], 164 Dy [7], and 168 Er [8] 
with large long-range anisotropic magnetic dipolar inter- 
action. Bosonic polar molecules with much larger elec- 
tric dipolar interaction are also being considered for BEC 
experiments [5j [9]. Thus one can study the properties 
of a dipolar BEC with variable short-range interaction 
[U [7] using a Feshbach resonance [10] . The dipolar BEC 
[TT] with anisotropic long-range atomic interaction has 
many distinct features pQ [5j [T2HT4] . The stability of a 
dipolar BEC depends on the scattering length as well as 
the trap geometry [4j[12j[14]. A disk-shaped trap, with 
the polarization z direction perpendicular to the plane of 
the disk, leads to a repulsive dipolar interaction making 
the dipolar BEC more stable [4]. On the other hand, a 
cigar-shaped dipolar BEC oriented along the polarization 
direction leads to an attractive dipolar interaction and 
hence may favor a collapse instability [4} H414T6] . Also, 
the anisotropic dipolar interaction leads, in general, to 
a distinct anisotropic density distribution in a dipolar 
BEC. The shock and sound waves also propagate with 
different velocities in different directions in a dipolar BEC 
[T7] . Anisotropic collapse has been observed and studied 
in a dipolar BEC of 52 Cr atoms [18]. 
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The properties of a BEC have been studied on different 
types of traps, such as, the harmonic trap pp, optical- 
lattice (OL) trap [19] , bichromatic OL trap [20], optical 
speckle potential trap [21] , double- well trap [22 , toroidal 
trap [23] [24] , ring-shaped trap in one [25] and three [26] 
dimensions, among others. To study the properties of the 
dipolar BEC, apart from the harmonic trap, the following 
traps have been used: OL trap [27 , bichromatic OL trap 
[28] , toroidal trap [29] and double- well trap [30] . 

Very recently, there has been experimental realization 
of ring-shaped and spherical-shell-shaped traps [3TJ [32] . 
The ring geometry was created by the time-averaged adi- 
abatic potential resulting from the application of an os- 
cillating magnetic bias field to a radio-frequency-dressed 
quadrupole trap [32] . The shell geometry was made from 
a cylindrically symmetric quadrupole field with its sym- 
metry axis aligned with gravity [32] . These geometries of 
the trap present an opportunity to study the superfluid 
properties of a BEC in a multiply connected geometry 
[32] not found in usual traps. In this paper, we study 
the properties of a dipolar BEC in ring- and shell-shaped 
traps. The shell-shaped trap is spherically symmetric, 
whereas the ring-shaped trap is at most axially symmet- 
ric. Different orientations of the ring-shaped dipolar BEC 
with respect to the polarization direction are considered. 
The stability of the dipolar BEC in these traps is illus- 
trated by phase plots of dipolar and contact interactions. 
The anisotropic nature of stability properties in a ring- 
shaped trap is further demonstrated by a consideration 
of the chemical potential of the system. The ring-shaped 
dipolar BEC is more stable when the plane of the ring is 
aligned perpendicular to the polarization direction, than 
when it is aligned parallel to the polarization direction. 
In the former case the dipolar interaction is mostly re- 
pulsive and in the latter case it is mostly attractive. The 
anisotropic density distribution of the dipolar BEC in 
these traps is explicitly demonstrated. 

We also consider a vortex dipolar BEC (rotating 
around the polarization z direction) of unit angular mo- 
mentum [33j [34] in ring- and shell-shaped traps, when 
the trapping potential is axially symmetric around the 
z direction. In the case of a ring-shaped trap the ax- 
ial symmetry is maintained when the ring is in a plane 
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perpendicular to the z direction. The ring-shaped vortex 
dipolar BEC is found to be nearly identical to the normal 
(nonrotating) dipolar BEC for the same sets of param- 
eters. The shell-shaped vortex dipolar BEC is found to 
possess a distinct density distribution when compared 
with a normal dipolar BEC. 

In Sec. II we describe the ring- and shell-shaped traps 
and present the mean-field Gross-Pitaevskii (GP) equa- 
tion which we use to study the normal and vortex dipolar 
BEC in these traps. In Sec. Ill we present the numer- 
ical results obtained by solving the GP equation using 
the Crank-Nicolson approach. We present stability phase 
plots of the dipolar BECs in terms of contact and dipolar 
interactions. The anisotropic density distribution of the 
BECs, which is a consequence of the anisotropic dipolar 
interaction, is illustrated for both the traps. The differ- 
ence in the density distribution between the normal and 
vortex BECs is also demonstrated. Finally in Sec. IV we 
present a brief summary and concluding remarks. 



II. ANALYTICAL CONSIDERATION 

A trap in the shape of a spherical shell will be taken 
in the form 



V(r) = V x -muj 2 (r -r ) 2 , 



(1) 



with r = {x,y,z},r = y/ (x 2 + y 2 + z 2 ), where m is the 
mass of an atom, ro is the radius of the shell, Vb is the 
strength of the trap and uj the frequency. For ro = 0, 
trap Q reduces to the usual harmonic trap. Similarly, 
the ring-shaped trap can be written as 

y(r) = y x^ W 2 [( v ^+?-ro) 2 + 9 2 ] ) (2) 

where p = (y cos a + z sin a) , q = (z cos a — y sin a) , Vo is 
the strength of the trap, a is the angle between the po- 
larization direction z and the perpendicular to the plane 
of the ring, and ro is the radius of the ring. For ro = 
and a = or 7r/2, trap ([5| reduces to the harmonic trap. 

A dilute dipolar BEC of N atoms, in shell- or ring- 
shaped traps (i]) or ([5]) will be studied using the following 
Gross-Pitaevskii (GP) equation [4} [T2] 



1 dt 



V + g\0\ 2 + 9dd F 



(M), (3) 



with g = AnaN^gdd — 3a^iV ". Here the dipolar non- 
linearity F = J U dd (r - r')\0(r' : t)\ W, U dd (H) = (1 - 
3 cos 2 0)/R 3 , R = r — r', normalization J (j)(r) 2 dr = 1, 
the angle between R and the polarization direction z. 
For the shell-trap 



V = V C 



shell 



^o(r-ro) 2 , 



and for the ring-trap 



V = V, 



ring 



V [(Vx 2 +p 2 -r ) 2 + q 2 



(4) 



(5) 



&dd — ^of^ 2 m/ (127rft 2 ) the strength of dipolar interaction, 
po the permeability of free space, and fit the (magnetic) 
dipole moment. In Eq. (|3|, length is measured in units 
of Iq = y/h/muj, potential and energies in units of hu, 
time t in units of to = . In this work we conveniently 
take Iq = 1 /im, which is the typical unit of length in 
BEC experiments. 

The ring-shaped trap with a = and the shell-shaped 
trap are axially symmetric around the z axis and a vor- 
tex BEC rotating around z axis with a conserved angu- 
lar momentum can be conveniently introduced in these 
cases. To obtain a quantized vortex of unit angular mo- 
mentum h around z axis, we introduce a phase (equal 
to the azimuthal angle) in the wave function [33]. This 
procedure introduces a centrifugal term l/[2(x 2 + 
the potential of the GP equation so that 



m 



V = V 



1 



2{x 2 +y 2 Y 



(6) 



where V is V^j^jj for the shell-shaped trap and ^ m g 
for ring-shaped trap with a = 0. For the ring-shaped 
trap with a ^ 0, the potential is not axially symmetric 
and a conserved angular momentum cannot be defined. 
We adopt this procedure to study an axially-symmetric 
vortex in a ring- and shell-shaped dipolar BEC. 

The chemical potential /i of a stationary state propa- 
gating as (j)(t) ~ exp(— ifit) is defined by 



/ 



; |V0| 2 + F|</>| 2 +^ 



■g dd F\^\ 2 



dr. (7) 



III. NUMERICAL CALCULATION 

We numerically solve the 3D GP equation ([3| using the 
split-step Crank-Nicolson method [35 j. The dipolar po- 
tential is divergent at short distances and hence the treat- 
ment of this potential requires some care. The integral 
over the dipolar potential is evaluated in Fourier (mo- 
mentum) space by a convolution identity [12] requiring 
the Fourier transformation of the dipolar potential and 
density. The Fourier transformation of the dipolar po- 
tential can be analytically evaluated [12] . The remaining 
Fourier transformations are evaluated numerically using 
a fast Fourier transformation algorithm. In the Crank- 
Nicolson discretization we used space step 0.1, time step 
0.002 and up to 256 space discretization points in each 
of the three Cartesian directions. We performed an anal- 
ysis of errors of sizes and energies calculated using our 
routine and find that the maximum numerical error in 
the calculation is less than 0.5 %. 



A. Normal dipolar BEC 

In Eqs. Q and ([5|, for shell and ring shapes, we take 
the strength of the potentials Vq = 10 and radius ro = 10. 
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FIG. 1: (Color online) Phase plot of contact and dipolar non- 
linearities g(= 4naN) and gdd{= SaddN) showing the stability 
line for shell- and ring-shaped dipolar BECs for a — and 
7r/2. The BEC is stable in the region below these lines and 
unstable above. 



FIG. 2: (Color online) The chemical potential \i versus the 
contact-interaction nonlinearity g for the ring- and shell- 
shaped BEC for different dipolar nonlinearity gdd- We show 
the results for gdd = and 200 for a = 0,7r/4, and tt/2. The 
arrow represents the onset of collapse in different cases. 



This makes reasonably strong traps, so that the widths 
of the shell or ring is small compared to the radius of the 
shell or ring and the shell and ring shapes of the BECs 
are pronounced. First we study the stability of the shell- 
and ring-shaped dipolar BEC. The stability properties 
of these BECs are best illustrated in phase plots involv- 
ing the strengths of contact and dipolar nonlinearities g 
and gdd as shown in Fig. [I] For the ring shape, we show 
two orientations of the ring for a = tt/2 and 0, the former 
corresponding to a ring in a plane parallel to the polariza- 
tion direction z and the latter in a plane perpendicular 
to the polarization direction z. The orientation of the 
dipoles in the former is quite similar to that in a cigar- 
shaped BEC along z axis in a harmonic trap, and that of 
the latter is similar to a disk-shaped BEC in x — y plane 
in a harmonic trap. The dipolar interaction contributes 
predominantly attractively in a cigar configuration with 
dipoles arranged parallel to the length of the cigar and it 
contributes predominantly repulsively in a disk configu- 
ration with dipoles arranged perpendicular to the plane 
of the disk. Hence for a ring-shaped dipolar BEC with 
a = the dipolar interaction contributes repulsively and 
to a positive quantity in energy and chemical potential. 
On the other hand, for a ring-shaped dipolar BEC with 
a = tt/2 the dipolar interaction contributes attractively 
and to a negative quantity in energy and chemical po- 
tential. Consequently, for a = 0, the ring-shaped dipolar 
BEC is stable for a reasonably large dipolar interaction 
below a critical value (gdd < 9dd^) without any contact 
interaction (g = 0). However, the dipolar BEC without 
any contact interaction (g = 0) collapses for larger dipo- 
lar interactions (gdd > 9dd^)- For a — tt/2, the dipolar 
interaction is mostly attractive and for a non-zero dipolar 
interaction, the ring-shaped dipolar BEC is stable only 
for contact interaction above a critical value as shown by 
the phase plot in Fig. [I] 



Next we consider the shell-shaped dipolar BEC. In this 
case the trapping potential is spherically symmetric and 
unless the dipolar interaction is large compared to the 
contact interaction, the contribution of the dipolar inter- 
action to the stability of the BEC, and to the chemical 
potential of the system, is insignificant. However, for a 
finite gdd the system collapses unless g is larger than a 
critical value. To avoid the collapse, a finite amount of 
contact repulsion is to be introduced through a finite g 
as shown in Fig. [I] 

Further consequence of stability characteristics is 
shown in Fig. |2j where we plot the chemical poten- 
tial of the dipolar BEC versus g for different values 
of 9dd(= 0,200). For a ring-shaped dipolar BEC with 
a = 7r/2, the dipolar interaction is attractive and con- 
tributes negatively to the chemical potential; whereas, 
for a = 0, the dipolar interaction is repulsive and con- 
tributes positively to the chemical potential. Conse- 
quently, the chemical potential for a = is the largest 
and for a = tt/2 is the smallest for a fixed g and gdd- The 
chemical potential for a = 7r/4 lies close to the nondipo- 
lar BEC with gdd = 0, as the contribution of the dipolar 
term is small in this case. For a shell-shaped dipolar 
BEC, as noted before, the contribution of the dipolar in- 
teraction to the chemical potential is insignificant and 
the chemical potentials for gdd = and 200 are prac- 
tically the same. The dipolar interaction energy nearly 
vanishes when integrated over all angles in a spherically 
symmetric density configuration [15] as in the spherical 
shell. However, in all cases of dipolar BEC (gdd — 200), 
the system collapses for g less than a critical value g c 
(indicated by an arrow in Fig. [2|, where the chemical 
potential suddenly jumps to a infinitely large negative 
value corresponding to a collapsed state. 

Next we illustrate the structure of a ring-shaped dipo- 
lar BEC by three-dimensional (3D) contour plots of den- 
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FIG. 3: (Color online) 3D contour plot of density \(p\ 2 of 
a ring-shaped dipolar 164 Dy BEC for add = 130cio,N = 
5000, a = 120a and (a) a = 0, (b)a = n/6, (c) a = tt/3, 
(d) a = 7r/2. The density at the contour is 0.005. The vari- 
ables x, y, and a are in units of Zo(= 1 /tm.) 



sity \<j>\ 2 . We exhibit the contour density plots of a 164 Dy 
BEC of 5000 164 Dy atoms with a dd = 130a 0. The S- 
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FIG. 4: (Color online) 2D contour plot of density in the x — z 
plane \<j>(x,y — 0,z)| 2 for a — 7r/2, for a ring-shaped dipolar 
BEC of 5000 52 Cr atoms with add — 15ao, and (a) a = 20ao, 
(b)a = 50ao, (c) a = 70ao, (d) a = 100ao. 




FIG. 5: (Color online) 3D contour plot of density |0| 2 for a 
shell-shaped dipolar BEC with gad — 100, and (a) g — 300, 
(b)g = 600, (c) g = 1000, (d) g = 2000. The density on the 
contour is 0.0005. The variables x, y, and z are in units of 
lo(= 1 /xm.) 



wave scattering length for contact interaction in this case 
is taken as a = 120ao in close agreement with some exper- 
imental estimate [7]. The 3D contour plots of density for 
a = 0,7r/6,7r/3, and ir/2 in this case are shown in Figs. [3] 
(a), (b), (c), and (d), respectively. For these parameters, 
in all the plots of Figs. |3j g = 399 and gad = 103. In 
the plane of the ring, the density is circularly symmetric 
for a — 0. However, as the angle a increases due to the 
anisotropic dipolar interaction the density distribution in 
the plane of the ring is no longer isotropic as can be seen 
in Figs. [3] (b) and (c). The anisotropy in density distri- 
bution in the ring is visible for a = 7r/3 in Fig. [3](c) and 
is most explicit for a = tt/2 in Fig. [3] (d) . Although, the 
dipolar interaction is long-range in nature, it is practi- 
cally zero across the diameter of the ring or the shell (20 
micron or 200,000 A). 

In the case of the ring the density is most asymmet- 
ric for a = 7r/2 and we study the density in this case in 
some detail for a dipolar BEC of 5000 52 Cr atoms with 
add = 15ao for different values of contact interaction. For 
this purpose, we show the 2D contour plot of density in 
the x—z plane y = 0, z)\ 2 in Figs. |4](a), (b), (c), and 
(d) for a = 20ao, 50ao, 70ao and 100ao, respectively. The 
anisotropic dipolar interaction is most pronounced when 
the contact interaction is the smallest, e.g., for a = 20ao 
in Fig. [I] (a). In this case the anisotropic distribution 
of density in the plane of the ring is most visible. As 
the isotropic contact interaction increases, the effect of 
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FIG. 6: (Color online) Densities |<y3(a;)| 2 = 
|^,0,0)| 2 ,|^)| 2 = |^(O,O,0)| 2 ,|^(r)| 2 = \4>(r,r,r)\ 2 , 
along x, z, and r radial directions for a shell- shaped dipolar 
BEC with g dd = 100 and g = (a) 300, (b) 600, (c) 1000, and 
(d) 2000. 



B. Vortex dipolar BEC 

Now we study the singly-quantized vortices in ring- 
and shell-shaped dipolar BECs. In such a vortex, the 
trapping potential should have azimuthal symmetry. In 
the ring shape, such symmetry exists only for a = 0. 
However, for a = the matter is localized far away from 
the center [viz. Fig. [3] (a)]. Consequently, the added 
centrifugal part in Eq. ^ is practically zero for large 
x and y 10) and the density of the vortex state is 
practically identical with the normal state shown in Fig. 
[3] (a). In the x — y plane, both have zero density at 
the center of the ring. Nevertheless, the phase of ring 
vortex wave function of small finite radius changes by 2ir 
as one moves in a closed contour around the center of the 
vortex. However, for larger values of angular momentum 
of the ring-shaped BEC there could be dissipation when 
the rotational velocity exceeds a critical value [32j [36] , 
larger than the velocity for the unit angular momentum 
considered here. Below this critical velocity the mass 
flow is dissipationless in analogy to electrical current in 
superconductors. Such a persistent flow was observed in 
a nondipolar BEC in an optically plugged magnetic trap 

Next we consider a shell-shaped dipolar BEC rotat- 
ing around the polarization z direction with unit angular 
momentum. For this purpose, we solve the GP equation 
with potential For the vortex state the density on 
the z axis |<^(z)f becomes zero. Otherwise, the added 



the anisotropic dipolar interaction is less and less pro- 
nounced and an almost symmetric density distribution 
in the plane of the ring is obtained for a = 100ao as can 
be seen in Fig. [4] (d). 

In Figs. [5] (a), (b), (c), and (d) we show the 3D con- 
tour plot of density |</>| 2 for gdd = 100 and different 
g = 300, 600, 1000 and 2000. For a clean visualization of 
the spherical-shell shape only one half of the full density 
distribution is shown. The density is most asymmetric in 
the polarization z direction in Fig. [5] (a) with the small- 
est contact nonlinearity g, while the dipolar interaction is 
most prominent. This asymmetry in density reduces as 
the contact nonlinearity increases making the anisotropic 
dipolar interaction less and less prominent. In Fig. [5](d), 
the density distribution is the most isotropic in 3D with 
dipole interaction playing a minor role for gdd = 100 and 
g = 2000. 

The above anisotropic density distribution of the shell- 
shaped dipolar BEC is further demonstrated by plotting 
the density along the x, z, and the radial r directions 
in Figs. [6] (a), (b), (c), and (d) for gdd = 100 and 
g = 300, 600, 1000 and 2000, respectively. The anisotropy 
in the density in the three directions is explicitly shown 
in these figures. The anisotropy clearly reduces as the 
contact nonlinearity g increases making the dipolar in- 
teraction less prominent, as can be seen in Fig. [6] (d) 
with the most isotropic density distribution. 
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FIG. 7: (Color online) Densities \(f(x)\ 2 = |0(x,O,O)| 2 and 
|c^(r)| 2 = |0(r, r, r)| 2 along x and r directions for a shell- 
shaped vortex dipolar BEC with gdd = 100 and g = (a) 300, 
(b) 600, (c) 1000, and (d) 2000. 
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centrifugal term in Eq. ([6| is very small for large x and 
y. For most parts of the shell the values of x and y are 
large and the centrifugal part of the potential has very lit- 
tle effect on the density of the condensate. Consequently, 
the density distribution of the dipolar vortex BEC will be 
very similar to that of the normal BECs shown in Figs. [6j 
with the only different that for the vortex the density at 
all points on the z axis will be zero: |^(^)| 2 = 0. This is 
illustrated in Figs. [7] where we plot the densities |^(x)| 2 
and |^(r)| 2 as in Figs. [6j It is seen that the densities in x 
and r directions of the normal and vortex BECs are quite 
the same. As the density distributions of the normal and 
vortex BECs are similar, the stability phase plot of Fig. 
[I] and the chemical potential plot of Fig. [2] for the shell- 
shaped normal and vortex dipolar BEC are practically 
the same. 

The difference in density distribution between 
the shell-shaped normal and the vortex BEC can 
be best studied by considering the isotropic two- 
dimensional (2D) density of the condensate \<fr(x,y)\ 2 = 
j dz\(j)(x, y, z)\ 2 obtained by integrating out the z depen- 
dence. For a vortex this 2D density is zero for x = y = 0. 
To illustrate the difference in density distribution of a 
normal and a vortex BEC, we plot in Fig. [8] the density 
|<l>(£, 0)| 2 versus x for the different cases shown in Figs. 
[6j For these sets of parameters the densities for the nor- 
mal and the vortex BECs are shown. The two densities 
are quite similar over most regions except near the cen- 
ter x = 0. For the vortex BEC the density is zero at the 
center, whereas for the normal BEC it has a small finite 
value. 



the stability plot in Fig. [T] 

The vortex dipolar BECs are also considered in the 
axially-symmetric ring- and shell-shaped traps. The ring- 
shaped trap is axially symmetric when the ring stays in 
a plane perpendicular to the polarization direction. In 
the case of the ring-shaped trap the density distribution 
of the normal and vortex dipolar BECs is practically the 
same. The difference in the density distribution of the 
normal and vortex dipolar BECs in the shell-shaped trap 
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IV. CONCLUSION 

We studied the properties of a normal (nonrotating) 
and vortex dipolar BEC in ring- and shell-shaped traps 
using the mean-field GP equation. The stability of the 
system with anisotropic density distribution, which is 
a consequence of the anisotropic dipolar interaction, is 
studied in phase plots of dipolar and contact interac- 
tions. The system is more stable with reduced anisotropic 
density distribution when the dipolar interaction is small 
compared to the contact interaction. In the shell-shaped 
trap the dipolar BEC is always unstable when the con- 
tact interaction is zero or attractive. In the ring-shaped 
trap the dipolar BEC can be stable for zero contact in- 
teraction when the plane of the ring makes an angle with 
the polarization direction (a ^ tt/2), as can be seen from 



FIG. 8: (Color online) Radial 2D density along x axis 
\<&(x, y = 0)| 2 = f dz\<j){x, y = 0, z)\ 2 of a shell-shaped normal 
and vortex dipolar BEC with gad = 100 and g — (a) 300, (b) 
600, (c) 1000, and (d) 2000. 



is carefully examined. The experimental realization of 
the shell- and ring-shaped traps [3TJ [32] and the present 
theoretical study will trigger further studies of dipolar 
BEC in these novel traps. 
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